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Abstract 

A new type of perturbation expansion in the mixing V of local- 
ized orbitals with a conduction-electron band in the U — > oo Anderson 
model is presented. It is built on Feynman diagrams obeying standard 
rules. The local correlations of the unperturbed system (the atomic 
limit) are included exactly, no auxiliary particles are introduced. As a 
test, an infinite-order ladder-type resummation is analytically treated 
in the Kondo regime, recovering the correct energy scale. An extension 
to the Anderson-lattice model is obtained via an effective-site approx- 
imation through a cumulant expansion in V on the lattice. Relation 
to treatments in infinite spatial dimensions are indicated. 
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The large Coulomb repulsion U on localized /- or (/-orbitals in strongly 
correlated electron systems is the source of exciting many-body effects and 
at the same time responsible for great technical difficulties: in a theoretical 
description based on Anderson- or Hubbard-type models, the most suitable 
starting point uses local (atomic) many-body states, which are perturbed 
by an intersite hopping t or a hybridisation V with delocalized band states. 
The atomic limit is apparently not given by a quadratic (i.e. one-particle) 
Hamiltonian H , and the use of a conventional Feynman-diagram technique 
in a perturbation expansion with respect to t or V poses an yet unsolved 
problem. 

In 'slave-boson' techniques Wick's theorem is recovered via auxiliary par- 
ticles, but for the price of a constraint [I, 2] to be maintained. It leads, if 
taken strictly [3] , to the unconventional time-ordered perturbation theory 
[4, 5] . Other approaches [6, 7] start from commutation relations of Hubbard 
operators [8, 5] and aim at an adequate 'Wick's theorem', which appears as 
a recurrance relation and does not provide with strict diagram rules [7] . 

We will show in the following for the U — > oo Anderson-impurity model 
[9] the possibility of an expansion in the hybridisation V using conventional 
Feynman diagrams. The strong local correlations are kept exactly with no 
need of 'slave' particles. As a test case for the new approach the Kondo 
problem is reconsidered, and finally an extension to the Anderson-lattice 
model is discussed. 

With Hubbard transfer operators [8] X aa = (1 — n^_ a )n^ a , 
X 0a = (1 - n S _ a )f ai X a0 = (X 0ff ) t , nl = flf a the U oo model is 
formulated as 

H = E £kA„ c k* + E 4x aa + E (v k x°°c kl , + Vkc{ a x°°) . (i) 

The effect of X-operators on Fock-space states \a; N c ) , characterized by 
an atomic /-configuration a £ {0,^,4,, 2} and some set A^ c of conduction- 
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band-occupation numbers, differs from that of canonical /-operators w a , f a , 
fl through the projecting prefactors: these allow for non- vanishing transi- 
tions only in the subspace without local double occupancy, X aa \(r'; N c ) = 
\<r;N c )8 a ,, a , X°°\<t';N c ) = |0; N c )8 a ,, a , X°°\0;N c ) = \a; N c ) , a =t,|. The 
trace in e.g. the partition function Z = Tr'[e~ f3H ] includes only this restricted 
Hilbert space (denoted by the prime) i.e. matrixelements like (M; N c \e~ pH \M; N c ) 
with M = 0, a . The X-operators create intermediate states according to the 
allowed transitions above. Since the trace offers only admissible local con- 
figurations the prefactors in X aa and X 0a are actually inoperative as long 
as the creation of a double occupancy is still prevented in keeping X a0 as it 
is: the replacement of X aa by and X 0a by f a does not alter the transi- 
tions occuring in the matrixelement. The term under the primed trace may 
therefore be simplified, leading to Z = Tr'[e~ f3H ] with 

H = H° + J2{v k X^c ka + Viclf a ) , H = T,ekA t , c k* + T,etflf* ■ 

In general this argument [10] holds for any matrixlement 

<M'; N' c \0{{X°°, X°% X"})\M; N c ) = <M'; N' c \0{{X°° , f a ,fj a })\M; N c ) 

(2) 

of any operator O in the subspace M, M' = 0,^,4, involved in the U — > oo- 
problem. Note, the matrix of a Hermiteian operator like O = e~^ H remains 
in this subspace Hermiteian. 

The local /-propagator F CT , CT /(r,r') = - < T{X° CT (r)X CT °(r')} > can now 
be written as F CT)CT /(r, r') = — Tr'[e~ f3H T{f a (T)X a °(t')}] / Z with (imaginary 
time) Heisenberg operators involving H according to (2). It still requires a 
reduced trace (Tr'), the quadratic unperturbed Hamiltonian H° alone is not 
sufficient in order to get a Feynman-type expansion instead of (time ordered) 
Goldstone diagrams [11]. But since X a °(t') = (/_ CT //I CT ;/J;)(r') eliminates 
all matrixelements (2; iV c | . . . |2; iV c ) , we may enlarge the trace to the full Fock 
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space and are left with a conventional interacting two-particle propagator: 

F,At,t') = ~ f <r{/4r)/-,(r' + + )/Ur')/i ( (r')}> 5 . (3) 

Wick's theorem is now applicable and an expansion in connected diagrams 
is formulated as usual (see, e.g. [12]) : a diagram of order V 2n consists of n 
simple hybridization vertices (Fig. 1 top) and n two-particle vertices (Fig. 1 
bottom). These are linked by bare conduction-electron lines G^ a {iui) = 
l/[iu>i — Ska] (full lines) and bare 'unprojected' /-propagators G°(r, r') = 
- <r{f a (r)fKr')} > So ^ G°(z^) = 1/M -ef] (dashed lines). Its 
sign ( — l) Uc is given by the number of fermion loops n c , and each internal 
frequency iu n is accompanied by a factor e - tuj ™°+ (the exponent differs in sign 
from standard rules !). 
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Figure 1: Vertices resulting from the application of Wick's theorem. Spin, 
frequency and direction of arrows are conserved. Dots represent matrix ele- 
ments Vk , V£ ■ The two-particle vertex (bottom) permits an f -creation only 
if a hole of opposite spin is present. 

The skeleton representation of the physical /-Green's function (3) shown 
in Fig. 2 involves fully renormalized 'unprojected' Green's functions G a {iui) 
(dashed double lines) and an irreducible vertex function T 2 [12] . A prefactor 
appears through 

< X 00 + X" > = < 1 - nL a >= ^Tr'[ e -^(1 - nL a )] = f < 1 - nL a > S , 
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it already reveals the exact spectral weight according to the sum rule 
IZo dupi(u) = < X 00 + X™ > , pf(u) = - fimF^iu + z0+) . The two 
fermion loops emerging from contracting three of the four external time labels 
of (3) contain the overall sign. In the one-particle part (Fig. 2 , hrst term) 
the loop also cancels the prefactor's denominator < 1 — nL a > H ■ 
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Figure 2: Skeleton representation of the physical f -Green's function (see 
text). Conservation of incoming frequency and spin is maintained by a 
dummy vertex. 

The Kondo problem is now considered within the perturbation scheme 
formulated above in order to identify important diagram types. The sim- 
plest reasonable result for the /-spectral function p^(oj) at low energies u is 
obtained by infinite-order resummation of spin-flip-scattering processes [13] , 
thus it may be expected that a simple mean-held (MF) theory is insufficient: 
based on a neglect of T 2 (Fig. 2 , second term) and the self-consistent treat- 
ment of G(iuJi) in lowest order (spin-degenerate) skeleton-self energy [14] only 
a slight shift of the bare /-level shows up. In going further, the analytical 
solution [15] of Bethe-Salpeter equations for the particle-hole (ph) channel 
and for the particle-particle (pp) channel, which contains the important spin- 
flip-ladder diagrams, produces via T 2 the Abrikosov-Suhl resonance at low 



5 



temperature T 



P' 




T-s-0 



ttT a 



A 



£(cu - Ta) , T A = Dexp 



( 



2A 



/ 



) 



including the proper energy scale T4 [13, 16] . The constant background due 
to the (ph)-channel has been left out here. D is half the width of a flat 
conduction band, A = vr | V | 2 / 2D at = V . 

Our result compares very well with analytical non-self-consistent resum- 
mation in time-ordered perturbation theory [17, 18] : in setting 
<X 00 + X aa >sa jj = \ the spectrum of [19] is exactly obtained. 

An extension of the method to the U — > 00 Anderson-lattice model, how- 
ever, remains difficult. The operator replacement (2) still holds for each 
lattice site = fi (as — )■ f ^ etc.), but in the local /-Green's func- 
tion F^^^iiui) double occupancies on all sites other than ji do contribute to 
an enlarged trace. If the replacement and trace-enlargement steps are per- 
formed after expanding F^^^iiui) in V first, a bookkeeping of lattice sites 
involved turns out to be necessary for each diagram, thus a linked-cluster ex- 
pansion required to get an intensive Green's function is prevented. The well 
known excluded- volume problem [5] showing up here can be treated by way 
of self-consistent resummations of lattice processes to an effective Anderson- 
impurity model [20, 21, 22] , a concept closely related to recent treatments 
of lattice models in infinite spatial dimensions [23, 24] . Rather than to aim 
at an effective-site approximation via the V-perturbation theory introduced 
here, a very general result is obtained following the lines of Ref. [20] : 

The full conduction-band Green's function G < j ta (iuJi) = 1 / [iui — ej ea . — 
Yife (iu>i)] of the Anderson-lattice model at = V and (via exact equations 
of motion) the T-matrix 
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(4) 
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are expressed through the band-self energy . The latter is expanded in V 
using bare cumulants [25, 26] which involve only the atomic (V = and any 
U) Hamiltonian and bare band-electron Green's functions . An effective- 
site approximation is made in the corresponding skeleton expansion [27] (see 
Fig. 3) through a decomposition of G% a = G c a + G% a , Q% = jj^; Efe G% a , 
Gfe = — Q c a and the neglect of its non-local part . 




Figure 3: Skeleton- cumulant expansion of the band-self energy for the Ander- 
son lattice. Full conduction-band lines G c kcr connect bare atomic cumulants 
(filled circles), which are k- (and a-) conserving and local, i.e. k -independent. 
They carry a at each incoming (outgoing) band line, also if amputated. 

Now S^ CT [Cr^/ ,] — y E^ mp is the local (fe-independent) band-self energy 
of an effective impurity model, and as a consequence the relations 

Gi = i/ [(Gir 1 - sr] , r ff = f ^E^ = s r / [i - Gi^r] (5) 

k 

hold by means of an expansion of Q c a and the local T-matrix T a in real space 
on the lattice. The propagator G° i a = emerging here is dehned as 

(7 flfl,<J' / j ,<7 l>\ l>\ ,<7 l>\ /i,<7 ' / j j!V\,7 V\ll\,<T V\V2,(7 l'2 lJ 2 ^ lJ 2l i ', <J 

with G° uu , a = ir 1 — Efe e lk ^- R ^G° h , E™? CT = 8„ ^E™ p and is according 

/i/i ,U # sites '—'n, 5 flfl ,<T A'jA 1 17 O 

to (5) to be identihed as the effective-impurity's 'bare' band propagator. It 
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can be determined from (4) and (5) resulting in a lattice T-matrix and self- 
consistency constraint as known from the XNCA ([22] , Eqn. (2.25), (2.21)): 

Tk, = 1 / [m- 1 - [GL - Gl\] , 1 = T7^ E 1 / [l " TAGL ~ Gl]} . 

(6) 

Since the expansion that leads to (6) containes no reference to any specific 
way bare (atomic) cumulants are to be calculated [26] , they were only used 
as a means, the 'outer XNCA-loop' (6) represents a general and well de- 
fined effective-site approximation. In this framework the present approach 
with 'unprojected' Green's functions in Feynman diagrams extends to the 
Anderson lattice. 

In the limit of infinite dimensions (d — > oo) [23] the local skeleton-self 
energy S™ p becomes exact following the arguments of Ref. [24] , which proves 
(6) an exact mapping scheme in this limit. Via equations of motion it can 
be shown to be equivalent to the 'usual' d — > oo-mapping scheme [23] and 
the one introduced in [28] . 

In conclusion, a way to set up an expansion of Green's functions around 
the atomic limit of the Anderson-impurity model via conventional Feynman 
diagrams has been demonstrated in considering the requirements for Wick's 
theorem without 'slave' particles introduced. Important diagram classes have 
been identified within an analytical treatment of the Kondo problem, indicat- 
ing the way to self-consistent conserving approximations. In order to obtain 
an extension to the Anderson-lattice model, however, we had to refer to the 
effective-site approximation by way of a general cumulant expansion around 
the atomic limit of the lattice. It remains the question whether correct spec- 
tral functions of the spin-degenerate Anderson lattice are thereby obtainable 
in the low temperature regime [22, 28] . This is under consideration via 
numerical investigation of conserving approximations mentioned above. 
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